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SUMMARY 

A numerical automation procedure has been developed to be used in conjunc- 
tion with an inverse hodograph method for the design of controlled diffusion 
blades. With this procedure a cascade of airfoils with a prescribed solidity, 
inlet Mach number, inlet air flow angle and air flow turning can be produced 
automatically. The trailing edge thickness of the airfoil, an important quan- 
tity in inverse methods, is also prescribed. 

The automation procedure consists of a multidimensional Newton iteration 
in which the objective design conditions are achieved by acting on the hodo- 
graph input parameters of the underlying inverse code. The method, although 
more general in scope, is applied in this paper to the design of axial flow 
compressor blade sections, and a wide range of examples is presented. 


INTRODUCTION 

Three dimensional flow analysis codes are increasingly playing a more 
important role in the design process of axial turbomachinery blading. These 
codes are generally used to analyze blade configurations that have been 
designed by two-dimensional methods. The fast computer generation of two- 
dimensional blade sections with prescribed aerodynamic characteristics still 
has a central role in the blade design process. 

Inverse hodograph codes for the design of controlled diffusion blades 
have proven to be an excellent source of innovative designs, both in the sub- 
sonic and the transonic regime (refs. 1 to 3). One advantage of these methods 
over direct design methods, aside from the generation of shock-free airfoils, 
is that the design philosophy can be incorporated into the design process from 
the start, by imposing a given surface speed distribution which reflects the 
desired aerodynamic behavior. Also, because the geometry is totally construc- 
ted by the method, it can produce body shapes that could hardly be predicted 
with direct methods (ref. 4). In general, the use of inverse hodograph meth- 
ods requires expertise on the part of the user; it hardly can be seen as a 
"black box" method. 

Both direct and inverse methods require an iteration procedure to gener- 
ate a blade section that meets the desired flow characteristics. In the 
direct design method a specified geometry is modified by interacting with a 
two-dimensional flow solver until the desired flow conditions are met. In the 
case of inverse hodograph methods, an airfoil with a prescribed speed distribu- 
tion is obtained but an iteration process is still necessary to achieve the 
other geometric and flow requirements. 


In this paper an automated procedure is presented that has been developed 
for use with the Inverse Hodograph Design Code (LINDES) described in refer- 
ence 5. In this code, three input hodograph design parameters and a prescribed 
surface pressure distribution control the flow characteristics at the design 
point. The solidity of the cascade, the inlet Mach number, and the inlet air 
flow angle are brought to their desired values by modifying the three hodo- 
graph design parameters in successive computer runs. The lift imposed by the 
input pressure distribution determines the air flow turning of the cascade. 

The automation procedure described in this paper consists of a multidimen- 
sional Newton method in which the objective functions are the cascade solid- 
ity, the inlet Mach number, the inlet air flow angle, the air flow turning, 
and the trailing edge closure or gap conditions. An input speed distribution 
with three free parameters is used in the automation procedure. The Newton 
iteration will then use six independent variables to achieve six objective 
functions . 


SYMBOLS 

dnte trailing edge thickness 
dste streamwise trailing edge gap 

M local Mach number 

Mq hodograph input parameter to control inlet Mach number 
M] inlet Mach number 

Qm peak value of input speed distribution 

Qte trailing edge speed 

R hodograph input parameter to control solidity 
Ste suction to pressure side arc length ratio 
inlet air angle 
A(J air flow turning 

0 hodograph input parameter to control inlet air angle 

£,n canonical complex characteristic hodograph variables 

p density 

a cascade solidity 

4>,Y potential and stream functions 
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FIGURES NOMENCLATURE 


BET1 

inlet air angle, [3] 


BET2 

exit air angle 


DELB 

air flow turning, A(3 


DNTE 

inviscid trailing edge thickness, dn^ e 


DSTE 

inviscid streamwise trailing edge gap, 

ds te 

LOSS 

loss coefficient 


Ml 

inlet Mach number, M] 


M2 

exi t Mach number 


SOL 

solidity, a 


TEV 

trailing edge thickness after boundary 

layer subtraction 

THCR 

maximum thickness to chord ratio 



THE BASE INVERSE HODOGRAPH METHOD 

The base inverse design code used for the automation procedure has been 
described in reference 1, and a users' manual was recently published (ref. 5). 
The purpose of this section is to review only those aspects of the basic 
method and design process needed to introduce the input parameters that are to 
be controlled by the automation procedure. A complete description of the base 
method is given in the above references and references therein. 

The inverse hodograph method constructs an analytical solution for the 
potential flow equations. The equations for the potential and stream func- 
tions cp, y, have the hodograph, canonical, complex characteristic form 
(refs. 1 and 2) 


^ = T + 


Cp = T U/ 

n - n 


( 1 ) 


where 


T+ = +1 


A 


7 


The computational complex variables £ and n are defined in a fixed, 
hodograph-1 ike, domain by means of a solution-dependent conformal transforma- 
tion. The details of the method have been described in references 1 and 2. 


The solution to the system (ref. 1) is expressed in the form 
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cp(£,n) = Re i cp 1 (^,n) log(n - n 1 > + <p 2 (£,n)log(n - n 2 ) + 

\ (2) 

t(£,n> = Re (5,n> log(n - np + 9 2 (£,n>log(n - n 2 > + 

The points n 1 and n 2 correspond to the location of the upstream and 

downstream singularities, respectively. They are located at the foci (ref. 1) 
of a fixed ellipse whose boundary corresponds to the subsonic part of the air- 
foil. The eccentricity of this ellipse determines the distance between these 
two points, and is prescribed in the code by means of a parameter R, the 
radius of a circular ring, 1 < |zj < R, topologically equivalent to the 
ellipse. By increasing or decreasing the parameter R, the solidity of the 
cascade a decreases or increases. The actual value of the solidity is part 
of the solution, hence the necessity of successive computer runs to adjust to 
the desired value. 

After the single valuedness of the stream function is imposed in 
equation (2) and the circulation over the airfoil is determined from the input 
speed distribution, two other arbitrary constants must be specified (ref. 1); 
this is done by means of two input parameters. These two parameters, Mq and 
©, control the inlet Mach number, M] and inlet air angle 3] and have to be 
modified appropriately until the desired values of Mj and 3] are achieved. 

It was stated in reference 1 that the elliptic transformation introduced 
there has the property that each of the input parameters R, Mq, and © have 
a dominant effect over each of the output parameters a, M] , and 3], respec- 
tively. It is this property of the elliptic mapping that is exploited in the 
Newton iteration process. 

As part of the solution, the body shape is found by using the integral 

/v e 

x + iy = J — - — (d 9 + i dy/p) . (3) 

The residue of this integral, given by the vector 

dx + i dy , (4) 

is another output quantity of the design code, it corresponds to the trailing 
edge opening. The automation procedure adjusts the normal and streamwise com- 
ponents of this vector until the desired trailing edge thickness dn^- e is 
obtained, while the streamwise trailing edge gap ds^e is forced to vanish. 
This produces a sharp trailing edge cut normal to the local flow direction and 
a prescribed trailing edge thickness. 


INPUT SPEED DISTRIBUTION 

During manual operation of the design code the design goals a, , and 
3] are achieved by acting on the input parameters R, Mq, and 0 previously 
described. Next, the input pressure distribution is adjusted to obtain the 
necessary flow turning A3 and the desired trailing edge gap condition 
(ref. 1). The flow turning is adjusted by varying the lift and the trailing 
edge gap by modifying the trailing edge speed and the relative arc length of 
the blade suction and pressure sides. 
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In order to be effective, the automation procedure must be able to use a 
family of pressure distributions defined by as few parameters as possible, but 
with sufficient versatility to produce the desired geometry. In this work, 
the speed distribution is defined by piecewise polynominal curves with all its 
coefficients except three fixed beforehand. For each type of application, dif- 
ferent polynomials can be chosen. In particular, Bezier polynomials have 
proven to be convenient for the present work. 

Figure 1 represents a generic speed distribution appropriate for axial 
compressor blades. In this curve, the trailing edge speed Qt e and the arc 
length Ste of either the suction or pressure side are left as free parame- 
ters on which the automated procedure can act. A third free parameter Q M is 
set as the difference between the peak speed on the suction and pressure sides 
and by controlling the area under the curve, it has a direct effect on the lift 
imposed by the speed distribution. 

It is worthwhile to note here that other parameters could also be left 
free, thus increasing the number of independent variables in the Newton iter- 
ation. For instance, the slope of the speed distribution at the leading edge 
has a direct effect on the maximum thickness of the airfoil, and could, con- 
ceivably, be added as a new free parameter. But, as stated earlier, the six 
free design parameters previously described have the important property that 
each has a dominant effect on one of the design objectives with only a second- 
ary effect on the others. 

The new form of definition of the speed distribution, which is essential 
to the automation procedure, has the additional advantage of making the design 
process less sensitive to spurious fluctuations in the speed distribution than 
when it was defined by a cubic spline, as it was previously done on the basic 
code. Besides, this new mode of defining the speed distribution does not 
seems to impose any limitation over the possibilities of the underlying code. 

In fact, because of the capability of making local changes in the areas 
needed, it promise to enhance these possibilities. 


AUTOMATION 

Consider the solution to the flow equations, generated by each cycle of 
the inverse code, as a vector valued function 

y = F<x> , (5) 


where 


x = (R.MQ.O.QM.Qte-Ste) (6) 

is the vector formed with N = 6 input parameters and 

y(o,Mi ,Bi .AP.dn^e .dste^ (7) 

is the vector of output conditions. The objective function for the automation 
procedure is the vector 
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y 0 = (a, Ml ,Pi ,A3,dn te ,ds te >0 


( 8 ) 


of specified output conditions. 

A Newton iteration is established to solve the vector equation 

F(x) - yo = 0 , (9) 


by means of the relation 


x n+l = x 

where J is the Jacobian matrix 


The numerical evaluation of the Jacobian matrix J requires N + 1 
cycles, one to calculate the starting point and N more, in which all but one 
of the independent variables are frozen, to compute the partial derivaties in 
J. Once the Jacobean matrix has been evaluated, it can be inverted and the 
new input vector x n+ ] is calculated. The process can then be repeated until 
a given tolerance 

Tol = max|(y - yo^il.i = 1»N (12) 


- J-l<x n Xy n - y 0 ) , (10) 



is achieved. 

Each iteration step consists of N + 1 cycles and yields a second order 
accurate method. We have found that faster convergence can be obtained by 
freezing the Jacobian J after the first N + 1 cycles have been performed. 
Although, in this case, the method is then first order accurate, only one 
cycle is required per iteration step. In either case, some relaxation is 
added for stability purposes. 

The method converges well, about forty cycles are enough to bring the tol- 
erance, 12, below a value of 10 -3 . In subsonic blade design the procedure is 
robust and, provided that a judicious choice is made, it will converge to the 
desired solution even from a relatively distant initial guess. 


EXAMPLES 

In this section, several examples are presented which have been developed 
to test the automation procedure. All the examples presented use the same 
input speed distribution shown in figure 1. As stated earlier, this speed dis- 
tribution is defined with only three free parameters, Q^ e , Ste» and Qm, on 
which the automation can act. The automation will search within this three- 
parameters family until a speed distribution is found that achieves the objec- 
tive function. In so doing, the speed ditribution shape will change according 
to the changes imposed on the parameters Qte. St e , and by the automated 
procedure. 
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One important point that the examples show is that starting from the same 
initial guess of speed distribution a broad range of airfoils can be obtained. 

A flat top speed distribution has been chosen, which is representative of a 
design philosophy. Accelerated profiles, for instance, could equally have 
been chosen and would result in different blade shapes. 

As a first example, a stator blade for a low speed fan is chosen. Four 
sections are designed, all with an inlet Mach number of 0.3 and inlet air 
angles of 45°, 40°, 35°, and 30° from hub to tip. The four sections, radially 
equidistant, have solidities varying from 2 at the hub to 1 at the tip. In 
all sections, the flow is requested to turn to the axial direction. 

The same initial guess is used for all sections, and an inviscid trailing 
edge thickness of 3 percent is imposed in all of them. Figures 2(a) to 5(a) 
show the inviscid airfoil and output surface Mach number distribution at the 
design point. Figures 2(b) to 5(b) present the cascade plane with the body 
shape after the boundary layer has been subtracted. 

The next example (figs. 6 to 9) consists of a higher speed blade, inlet 
Mach number of 0.5, and with the same air inlet angles, flow turning and solid- 
ities as the previous case. All the objectives were achieved, but the hub sec- 
tion presented an excessive amount of diffusion when a trailing edge thickness 
of 3 percent was specified. This section was then redesigned, imposing an 
inviscid trailing edge thickness of 3.5 percent. 

Finally, a transonic, shock-free, tip section is designed with the same 
flow conditions as the tip sections in the two previous cases, but with an 
inlet Mach number of 0.7 (fig. 10). Although the design is achieved, it 
required an initial guess closer to the design point. The parameter Mq has 
to be chosen well under the value that produces the objective design, in order 
to obtain first a subsonic blade of approximately the same solidity, inlet air 
angle and flow turning. The shock-free design problem is intrinsically differ- 
ent from the subsonic design and special care has to be exercised to produce a 
relevant case. Although the method solves a mathematically well posed prob- 
lem, see reference 1 and references therein, there is no guarantee that every 
given pressure distribtuion will produce a physically meaningful shock-free 
solution, hence the necessity of starting from a reasonable guess. 


CONCLUSION 

In this paper, the automated design of controlled diffusion blades is 
proven to be feasible. For this purpose, an automation procedure has been 
developed for the design of axial turbomachinery blading to be used in conjunc- 
tion with an inverse hodograph design code. With this automated design method, 
the inlet Mach number, inlet air flow angle, air flow turning, solidity and 
trailing edge thickness of a cascade of airfoils are prescribed. 

Applications to axial compressor stator blade design are presented. Two 
subsonic blades, each composed of four sections, and one transonic blade sec- 
tion are presented. All cases are produced with a single input speed distribu- 
tion, defined by three free parameters, on which the automation procedure acts 
to yield the different blade sections. 
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FIGURE 1. - GENERIC SURFACE SPEED DISTRIBUTION. 
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HI 0.299 
BET 1 29.99 
H2 0.263 
BET2 0.03 
DELB 29.97 
SOL 1.00 
THCR 0.099 
ONTE 0.030 
OSTE -0.001 
TEV 0.020 
LOSS 0.018 

(a) ( 
FIGURE 2. - LOW SPEED STATOR TIP SECTION. 






HI 

0.299 

BET 1 

35.01 

M2 

0.260 

8ET2 

0.09 

OELB 

39.97 

SOL 

1.19 

T HCfi 

0.112 

ONTE 

0.028 

OSTE 

-0.001 

TEV 

0.018 

LOSS 

0.02S 




(a) 


FIGURE 3. - LOW SPEED STATOR MID-TIP SECTION 
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<a> <b) 

FIGURE 4. - LOW SPEED STATOR MID-HUB SECTION 
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(a) 

FIGURE 5. - LOW SPEED STATOR HUB SECTION. 
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(a) (b 

FIGURE 6. - MID SPEED STATOR TIP SECTION. 
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(a) 


FIGURE 7. - MID SPEED STATOR MID-TIP SECTION 
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TEV 0.018 
LOSS 0.0*0 
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MID SPEED STATOR MID-HUB SECTION 





Ml 

0 . *9 8 

BET 1 

*5.26 

M2 

0.3*9 
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(a) 

FIGURE 9. - MID SPEED STATOR HUB SECTION. 
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(a) (b) 

FIGURE 10. - TRANSONIC STATOR SECTION. 
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